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Recently, quantum states of ultra-cold neutrons in the Earth's gravitational field have been ob- 
served for the first time. From the fact that they are consistent with Newtonian gravity on the 
10%-level, analytical limits on a and A of short-range Yukawa-like additional interactions are de- 
rived between A = 1 /ira and 1 mm. We arrive for A > 10 /im at a < 2 ■ 10 11 at 90% confidence 
level. This translates into a limit g s g p /hc < 2 ■ 10~ 15 on the pseudo-scalar coupling of axions in the 
previously experimentally unaccessible astrophysical axion window. 
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I. INTRODUCTION 



X 



A gravitational bound quantum system has been re- 
alized experimentally. In this experiment, ultra-cold neu- 
trons (UCN) are confined in between a bottom mirror and 
the gravitational potential of the Earth [3, H, Hj] . The neu- 
trons are slow enough that they are reflected from the mir- 
■ror at all angles of incidence. Therefore, the mirror can be 
modeled by an infinite high potential step. The neutrons 
are found in discrete quantum states of the gravity poten- 
tial. Between the UCN source and a UCN detector one 
places a quantum state absorber at a certain hight above 
the mirror. No neutrons except those in sufficiently low 
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FIG. 1: 

Experimental setup: Neutrons pass through the mirror- 
absorber system. The Airy functions for the first two bound 
quantum states are shown. The absorber is at variable height. 
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quantum states as given by the absorber height can pass 
through the slit between the mirror and absorber, and 
higher, unwanted states are removed and scattered out of 
the experiment (see Fig. [TJ. 

A side-effect of this experiment is its sensitivity to 
additional short ranged forces at length scales below 
10 jum 4, 5], while all electromagnetic effects are extremely 
suppressed compared to gravity The quantum states 
probe Newtonian gravity between 10 -9 and 10~ 3 m and 
the experiment places limits for gravity-like forces there. 
So far, significant limits on hypothetical forces from this 
experiment are mainly obtained from measurements of the 
ground state or from a fit to the data [3, [j| . Other experi- 
mental limits on extra forces are derived from mechanical 
experiments and can be found, e.g., in [1, 0,11, SIS- In the 
light of recent theoretical developments in higher dimen- 
sional field theory 0, [II [II [IJ 01 , gauge fields could 
mediate forces that are 10 6 to 10 12 times stronger than 
gravity at sub-millimeter distances (depending on the size 
of the extra dimensions), exactly in the interesting range 
of this experiment and might give a signal in an improved 
setup. Recent theoretical developments support this orig- 
inal idea of strong forces with bulk gauge fields. Burgess 
et al. [ill predict deviations from Newton's law on the 
micron scale on the basis of supersymmetric large extra 
dimensions (SLED). The basic idea behind this proposal 
is to modify gravity at small distances in such a way as 
to explain the smallness of the observed cosmological con- 
stant. A radius R of 10 microns as well as the necessary 
interaction strength may turn out to be well-motivated. 
We presents limits for additional interactions using an el- 
ementary particle, the neutron. The limits are derived in 
regions of the parameter space, where these interactions 
can be treated perturbatively. 
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II. ADDITIONAL SHORT RANGE FORCES 



We begin with recalling the standard parametrization 
of a fifth force by means of a Yukawa potential describing 
the low energy limit of the exchange of massive particles, 
proceeding then to the action of a perturbatively weak 
Yukawa-like fifth force on the gravitational bound states 
in different regimes of its range A. On the assumption that 
the form of the non-Newtonian potential is given by the 
Yukawa expression, for a source mass m and distance r 
the modified Newtonian potential 4>{r) has the form 



(r) = -G 4 — (1 + a-e 
r 



(1) 



where A is the range over which the corresponding force 
acts and a is the strength normalized relative to Newto- 
nian gravity. G4 is the gravitational constant. The mass 
of this extended source will modify the Newtonian poten- 
tial if strong non-Newtonian forces are present, and this 
can be seen if neutrons are present. For small distances 
z from the mirror, say several micrometers, we consider 
the mirror as an infinite half-space with mass density p. 
By replacing the source mass m by J pdV , the Yukawa- 
modification of the potential <fi(r) has the form 



A<f>(z, A) = -2tt • pa\ 2 Gi ■ e 



-z/X 



(2) 



Thus, the effective gravitational potential close above 
the mirror (close means: heights z « D, D: diameter of 
the mirror) is given as: 



9-z 

<l>0 



-2ir ■ a ■ A ■ G4 ■ p ■ e 



-z/X 



(3) 



where p denotes the mass density of the mirror material 
(glass in our experiment). 

The absorber has (besides its Fermi pseudopotential) 
in presence of a fifth force a Yukawa-like additional poten- 
tial attached to it just like the bottom mirror. Then A0 
becomes 



A(f>(z) = -2ir-a-\ 2 -G 4 - ( Pl ■ e' z / x + p 2 ■ e -(*-*)M) (4) 

where h denotes as before the absorber height and p\ , p 2 
the mass densities of the bottom mirror and the absorber, 
respectively. In some parts of this paper, we can ignore 
the absorber. 



Within an unchanged Newtonian gravitational poten- 



tial 



(j> {z) = g- z 



(5) 



the bound eigenstates of the UCN are given by Airy func- 
tions Ai(z). This function behaves similar to a sine wave 

below and approaches zero exponentially above this 
classical turning point height. With 

2 I 4 



-1 2/3 



(o) 



R = 



/ fi 2 \ 1/3 / 

( 9V J . C n = (nf 
\2migJ \ J 



-1/2 



one finds 



,(0)\ 



"9 I 



dz-^}Az)\z) , (z' \z) = 2tt • S(z' - z) 



V^M = C n -Ai((z-hW)/R),C n 



1? 



Cn- d -§{-h^/R)-^+0{{z/Rf) 



Cj7 



l/^\ 1/4 * 



R 



R 



■0((z/Rf) 



(7) 



(8) 



Here, denotes the position of the last turning point of 
the bound state ip^l(z), which coincides with the turning 
point height of the classical motion of a particle with en- 
ergy = mghn^ ■ hn^ has been determined [j| for the 
first two states to be 



, cxp 



: 12.2±0.7 stat ±l. 
21.6 ±0.7 sta t± 2.2 



'syst 
syst • 



(9) 
(10) 



With this analogy in mind, is considered as the 
height of the wave function. The ground state is opaque 
for neutrons as long as it sufficiently overlaps with the 
absorber above the mirror. This is not the case when the 
absorber position exceeds the height h„ following the ex- 
ponentially decaying tail of the bound state's Airy func- 
tion, and neutrons in that state are transmitted to the 
detector, see Fig. 1. In the WKB approximation the ipn}g 
are: 




(11) 



III. OUR METHOD 



We start from the observation that at macroscopic 
absorber heights of several microns, where light is easily 
transmitted, no neutrons pass through the gap between 
the mirror and the absorber [3[. A neutron in state n is 
only transmitted as long the absorber height is higher than 
h {0] 

1 in ■ 

Taking additional gravity-like forces into account, 
first order perturbation theory predicts a shift of the en- 
ergy eigenvalue of the ground state. The effect of the 
Yukawa term in eq. §5§ in first order is 



AM 1 



,(0) 



A<p(z) 



4 0) > 



(12) 



This energy shift is accompanied by a shift of the turning 
point by some additional height Sh, in this way chang- 
ing the onset of the transmission in the experiment. Our 
measurement follows the Newtonian expectation. 
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The absorber consists of glass material, where the sur- 
face has an approximate gaussian roughness of 0.75 /im. 
The height of the absorber has been calibrated with wire- 
spacers of known thickness, a mechanical comparator, and 
a long focus microscope. An absolute height calibration of 
better than 0.5 fj,m has been achieved. Next, we recall that 
the absorber forms a hard wall for the very slow neutrons. 
This results in a squeezing of the bound states [Tt], H3| 
and a small energy shift, which is considered in the follow- 
ing way: The squeezing increases the bound state energy 
eigenvalue by an amount Ai? squocz — m ■ g ■ Ag q , where 
As q « 0.8 [im denotes the corresponding shift of the turn- 
ing point of the squeezed wave function. As a result, the 
turning point of the squeezed bound states in presence of 
a real absorber is given by h n — hn + Ag q . We now calcu- 
late the influence of a hypothetical fifth force on the turn- 
ing points of the first two states (Eq[S]) and compare them 
with the experimental results given in Eq[51 The actual 
neutron transmission of the n bound state is observed 
at a mean absorber height h n + At ■ At has been chosen 
to be 3/Ltm, since above the turning point, the slope 
of the transmission curve shows that state n is transmit- 
ted. Summarizing both effects, a bound state n transmits 
neutrons, which are clearly visible in the detector if the 



implies a turning point condition given by 



absorber is at a height h n + Ax = h n 



(o) 



In this paper, we derive limits on 5th forces, which 
are largely independent from a precise knowledge of the 
absorber height and unaffected by the features of the ab- 
sorber. The limits are obtained from the fact that the 
ground state and the first excited state wave functions 
are differently affected by non-Newtonian forces. The fifth 
force would introduce the height shift Sh, so that we mea- 
sure h = hn^ +A.T + As q + Sh. In particular, the difference 
of the turning point heights of the ground state and the 
first excited state 



Ah^ :=h { °>-h { °> «10.3 M m 



(13) 



is consistent with the measurement Ah cxp = 9.4± 1.2 fim. 
This error includes in addition to the statistical uncer- 
tainty a 0.5/iin calibration uncertainty [|[ and a 30% un- 
certainty due to model dependence. (N.B.: Only the rela- 
tive error of a height measurement enters.) Since this rel- 
ative quantity is mostly insensitive to absolute offsets of 
the measurement process of h, we will use it to derive lim- 
its on the presence of additional short-range fifth forces, 
i.e., we can exclude a fifth force induced shift 

\5Ah\ > 1.64- cr Ah = 2.0 /xm (14) 

at 90% confidence level. In the case of Eq. (j2"9"|) , the choice 
of At and Ag q enters, but our limits are largely indepen- 
dent from these quantities. 

To apply this method it is essential to derive the con- 
nection between the perturbative correction to the energy 
eigenvalue of a given bound state, which is induced by the 
additional force, and a possible accompanying shift 5h n 
of the turning point of the n th state's wave function. The 
stationary Schrodinger equation for the bound states 

d 2 2m 

-Q-2^n,g{z) = --p- -[En-m- <j){z)} 1pn,g{z), (15) 



dz 2 



= 

= ■ [E n - m ■ 4>(h n )]ip n ^(h n ) 

=>■ E n — m ■ 4>{z) = for z — h n . (16) 



Write now for the height of the wave function under the 
presence of the additional force h n ~ h n + Sh n + As q , 
E n = E n 0) +AE^ ] +AE SQUCCZ and use the fact that E n 0) = 



nig ■ hffl ■ Then one arrives at 



AM 1 ) 



h=h n +Aq 



Sh n + A0(/4 O) + Sh 







A Sq ) 

J (17) 

with A<j>(z) given by eq. ([3]). Note, that the piece linear 
in As q from the squeezing of the bound states between 
absorber and mirror cancels against the corresponding en- 
ergy shift Ai? squccz . It is immediately clear from here that 
concerning the situation with just a bottom mirror for 
A << one has A^hn^) w and eq. (fT7|) simplifies in 
this regime to 



Sh„ = 



,(0) 



A<j>{z) 



,(o) 



(18) 



In the opposite case A > hn one may linearize Acf>(z) in 
z/A, which approximation then has to be used simulta- 
neously in exploiting eq. (fT7|) and eq. (jl"2]) to calculate 
AE^. 

For later use, let us note one further property of the 
above turning point condition. A<j)(z) may contain a con- 
stant, position-independent part (for instance, the above 
linear approximation in the case A j> generically pro- 
duces such a constant piece). Now, from considering the 
general behaviour of the exact solution it is clear, that 
changing the potential by an arbitrary constant must leave 
the whole bound state as well as its turning point un- 
changed though it does change the energy eigenvalue of 
the state. This fact is clearly contained in eq. (fT7|) : Imag- 
ine adding a constant A(f> const to A<f>(z). Then its contri- 
bution to AEn^ is given by A(j) const and thus cancels out 
against the same term in A(f)(z). Thus one may expand 
A(j)(z) around any given convenient point z and drop the 
constant piece. 

In either case we proceed then by extracting the turn- 
ing point shifts (5/io,i from eq. (jTTJ) for the ground state and 
the first excited state, respectively. Forming the difference 



5 Ah = 5hg — Shi 



(19) 



allows us then to extract limits on the strength a of the 
additional fifth force as a function of its range A by de- 
manding eq. ([TJ}, the experimental constraint at 90% C.L. 
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IV. BOTTOM MIRROR AND NO ABSORBER - 
CASE I: SMALL A < h { ° } 

Consider now the first case A << h^K Here, for a 
perfect absorber as said above, eq. (|T8|) provides a good 
description of the shift of the ground state turning point 
(see the general discussion of Sect. III). Then the linear 
approximation of the bound states given in Sect. II suffices 
to calculate eq. ([T5]) since A(f> then is confined to a region 
,(o) 



the non-penetration region and thus that feature of the 
measured neutron transmission function which responds 
most sensitively to energy or equivalently /i„-shifts, we 
evaluate eq. (fT5)l for n = 0: 



Sh n 



1 



,(0) 



A0 



G„ 



R 



a ■ A 2 • G 4 ■ p 
2-q-R 3 



dz ■ z 2 e- z ' x 



-Cl 



h (0) 

'hi 

~R 



A 5 ■ G 



4 • P 



g-R? 



(20) 



Forming SAh according to eq. (HHJ) we demand the ex- 
perimental constraint eq. Q14p . Therefore we arrive at an 
exclusion limit in a-A-space given by 



a < 



9-R 3 



\r<2. h-o n 2 \/B. i 



tt • A 5 • G 4 • p 



6 Ah ~ A" 



(21) 

which behaves symmetrical for attractive and repulsive 
forces. 



BOTTOM MIRROR AND NO ABSORBER 



CASE II: LARGE A >> h. 



(0) 



In the second case A >> Kq . Then A<j> ~ const. 

over the whole range where ipn} g is sizable. However, since 
constant pieces of the potential drop out from the turning 
point condition one has to apply it now in its precise form 
carefully expanding the fifth force potential to linear order 
in z and ho, respectively. This yields 

A<j>(z) = 2tt • a • A 2 • G 4 • p ■ (l - J) + O {z 2 /X 2 ) . 

Instead of A<f> we plug its contribution linear in z into 
eq. (fT2")) (recall that constant pieces of A<j> later will drop 
out of the turning point condition anyway) 



AEP = 



,(()) 



A<j> 



-Cl ■ m ■ 2tt • a ■ A 2 • G 4 • p 



dz ■ 



m ■ 2ir ■ a ■ A • G 4 • p ■ (z) 7 



(22) 



where eq. (jTTJ) allows to compute the turning point shift 
Sh of the ground state given by 

2n ■ a ■ A • G 4 • p (<*}„ - (^°) + A Sq )) 

- (g + 2?r • a ■ A • G 4 • p) ■ 5h n = . 

From this follows - via forming the quantity 8 Ah again - 
a limit in the case of a repulsive interaction (a < 0) given 

by 



A << h . Since it is the ground state which defines a > 



1 



2tt • G 4 • p i . Ah(o>-«s) -<*h) A 

r 1 ^ SAh 



A" 1 . (23) 



For the attractive case (a > 0) a smooth solution of 
eq. (HZJ) exists for all > SAh > -(Ah^ - ((z)o - 
which yields a limit 



a < 



1 



2tt ■ G 4 ■ p Ah(°)-(( z )o-{ z ) 1 ) 



- 1 



--A- 1 . (24) 



A 



Here it is (z) w 1.56 • i? w 9.15 pm and (z) 4 w 2.73 ■ R w 
16.0 /iTO. 

The transition between the two regimes of small and 
large A takes place just around A ~ 5 ... 7 pLtn as it can be 
seen by comparing eq.s (j2T|) and (j23j) . 



VI. BOTTOM MIRROR AND REAL ABSORBER 
- LARGE A 



A real absorber now has (besides its Fermi pseudopo- 
tential) in presence of a fifth force a Yukawa-like additional 
potential attached to it just like the bottom mirror. Then 
A(j> becomes 

A<f>(z) = -2tt • a ■ A 2 • G 4 ■ ( Pl ■ e- z ' x + p 2 ■ e"^-*)/^ 

where h denotes as before the absorber height and pi , p2 
the mass densities of the bottom mirror and the absorber, 
respectively. For large A it makes sense to expand eq. (|23|) 
in (z — h/2)/X around z = h/2: 

A(j){z) ~ -2vr • a ■ A 2 • G 4 • e"' l/2A • ( Pl + p 2 

z-h/2 ( Pl + p 2 ) ■ (z - h/2) 2 



(Pl - Pi) 



2 ■ A 2 



= const - tt a Gi ■ e~ h/2X ■ [ - 2A • (pi - p 2 ) ■ z 
+ ( Pl +p 2 )-(z 2 -z-h)]+X 2 0(^) . (26) 



A. Special case pi — pi 

Now in our case both mirror and absorber are made 
from glass, so p\ = p 2 = p. Then the direct linear term in 
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the former expansion vanishes - and thus all terms contain- 
ing positive powers of A! For large A thus just the quadratic 
term remains since it is independent of A, all higher terms 
are suppressed by negative powers of A and vanish in this 
limit. As an immediate consequence this means that for 
large A the presence of an absorber with a density equal 
to that of the mirror leads to a A-independent limit on the 
strength a of an additional Yukawa-like interaction. 

To determine this limit one computes the energy cor- 
rection eq. (fT3i|) induced by the potential eq. (p?fi|) now to 
first order in perturbation theory 

AEW = -m ■ 2na ■ G 4 ■ p ■ e~ h/2X ■ (z 2 - z ■ h) n , (27) 



where (z 2 - zh) n = (z 2 - zh ( n >) n - (z) n (A T + A Sq + 5h n ) 

for h — hn + oh n + At + Ag q denotes the expectation 
value of the z-dependent part of the correction (for the 
z-independent part dropping out of the turning point con- 



dition eq. (JTTJ) see the general discussion above) to the 
potential with respect to the ground state V'og- Then one 
inserts eq. (|26l) into the turning point condition eq. (|17p 
carefully evaluating to linear order in 5h n , uses the value 
of the turning point h n — h„ + 5h n + Ag q in the experi- 
ment, and obtains 



AEi 1 



h=h n +Ai 



— m ■ g ■ Sh n 



in 



M 0) + Skr. 



^Sq 



(28) 
. 



Determining 5h n from these equations for n = 0, 1 allows 
us to form subsequently SAh = Sh — Shi. Demanding 
then — 1.64a Ah < SAh < 1.64a Ah yields then a limit on 
the interaction strength a given by 



J 



\a\ < g\SAh\ 



it ( e -& [(z 2 ) + A T (4 0) + A Sq ) - <^)o(C + At + A Sq ) 



! + A T (hf } + A Sq ) - {z)x{hT + A T + A Sq 



(0) 



+5 Ah (e-^(A T - (z) ) + e _ ^(A T - {z)ij) +Vd) pG A 



(29) 



D = 



{e-^ [(z 2 ) Q + A T (C + A Sq ) - <z)o(4 0) + At + A Sq )] - [<z 2 )i + A T (/ l < 0) + A Sq ) 
-(^)i(4 03 + A T + A Sq )l + e-^(A T - (z) )SAh + e -^(A T - (z^SAhY 



-4e 



fiQ + hi f 

Ah i (Ai 



(z 2 ) + A T (4 0) + A Sq ) - (^)o(4 0) + At + A Sq ) 
+ (A T - (z) )5Ah] - (A T - (z) ) [(z 2 ) 1 + A^h^ + A Sq ) - (z)^ + A T + A Sq )] } 

I 



which is plotted in Fig. [5] These limits become essentially 
independent of A for A > h w hf\ The result for the 
repulsive case contains a pole, which in the worst case, 
e.g. for A = 500 pm is at 



SAh w 2.2 pm 



(30) 



and signals a breakdown of the perturbative approach. 
The pole in the attractive case above corresponds to the 
fact that for an attractive Yukawa potential at the mirror 
there is a \Sh\, above which the wave function begins to 
get sucked into the Yukawa potential. However, once the 
attractive potential gets strong enough for this to hap- 
pen, the perturbative expansion around the original states 
ceases to be a good approximation, which explains such a 
pole: If one added the higher orders of perturbation theory 
to the above constraint, the above pole would occur at just 
that Sh, where the corresponding wave function becomes 



non-perturbatively deformed, which would then coincide 
with validity boundary of perturbation theory derived be- 
low in eq. 



Finally, for really large A the constant term in the 
potential eq. (|26|) eventually becomes energetically domi- 
nant. Since this term in any case has to stay significantly 
smaller than the total kinetic energy Ejjc n of the ultra- 
cold neutrons (velocity of ~ 5 m/s) entering the waveguide 
(otherwise they cease to be transmitted entirely, which is 
not observed), the condition 



A(j)(z) ~ -2tt • a ■ A 2 • G 4 ■ p < 0.1 • E, 



UCN 



limits in the case of an attractive interaction the validity 
of eq. (j29f to A < 10 mm. For the case of a repulsive in- 
teraction this leads to a further decrease of the bound of 



G 



the strength a above A ~ 10 mm of the 5th force like 



\a\ < 3.8- 10 12 



10 mm \ ' 



(31) 



B. 



General case pi 7^ p2 



Let us now discuss what happens if pi 7^ pi. In this 
case the linear term in eq. (|26p reappears. Since this term 
is ~ A it will dominate the quadratic term considered so 
far for large A. Thus, if p\—pi is not tuned to be very small, 
the limit on the Yukawa interaction is again given through 
calculating 8 Ah by the eq.s |23|) and (j24|) . however, with 
p replaced by p\ — P2 and A<p given by the piece linear in 
z of eq. (HU 



< 



1 



2tt ■ G 4 • (pi - pi) Ah0)-« Z ) -( Z ) 1 ) , A 

SAh 



(32) 



A - , A > h 



Let us shortly note now, that this limit for large 
A > 20 pm can be relatively easily converted into bounds 
of the strength of the matter couplings of axions. Axion in- 
teractions with a range within 20 pm < A < 200 mm (cor- 
responding to axion masses 10 -6 eV < m a < 10~ 2 eV), the 
'axion window', are still allowed by the otherwise stringent 
constraints posed by cosmological data (see e.g. [2ll |22||). 
They lead to a potential which is proportional to the 5th 
force potential eq. ([2]), however, the axion- induced poten- 
tial changes sign with the direction of the neutron spin po- 
larization relative to the mirror. Thus, the relevant limit 
is again given by eq. ([55| with, however, one small but sig- 
nificant change induced by this pseudoscalar nature of the 
axion interaction, namely, that pi — ► — pi (see [23J]). This, 
in turn, implies that the bound on the scalar-pseudoscalar 
axion interaction could be derived using just eq. (|3"5|) , how- 
ever, with P2 replaced as p-x — > — p-x (or the perturbative 
limit eq. (|32p , with again the replacement P2 — * —pi)- 



P2 — > —pi - axion limits in the astrophysical axion 
window 



Here we have defined h = (hi + Hq)/2 ps 18.9 p and ne- 
glected factors exp(±A/i/2A) which are suppressed for 
A > Ah = —Aha/2 sw 5.2 p. This treatment is pertur- 
bative. 

We can cross-check this result by looking again at the 
potential eq. ([21?]). and add its linear piece to the earth's 
gravitational potential cf>o — g ■ z 

4>(z) = g ■ z + 2ir ■ a ■ X ■ G4 ■ (pi — p^) • z . 

P2 amounts to a renormalization of 



Clearly, the case pi 
g given by 



9 



9 + 8g 

Sg = 2ir ■ a ■ A • G4 • (pi — P2) 



(33) 



Since the turning points of the bound state Airy func- 
tiosn are given by eq. ([7]) the turning point shifts in this 
situation, where g renormalizes, using eq. (|33[) write as 



Sh n = 



3tt 
~2 

3tt 
T 

-h^ ■ 



2/3 



2/3 



dR . 
-OS-' 69 

R *9 



2n ■ a ■ A • G 4 • (pi — p 2 ) 
~9 



(34) 



This relation then implies after forming 8 Ah = 8ho — Shi 
through the experimental constraint eq. (fl4|) \8Ah\ < 
1.6 pm a symmetrical limit on the strength of the Yukawa 
interaction 



a < 



3 -9 



\8Ah\ 1 
2tt • G 4 • (pi - P2) ' AhM ' A 



(35) 



which is the same functional dependence as given in 
eq. (Aft(°> - ({z) - (z) n ) S3 A/i(°) /3) for large \>h. 



To make the last statement more precise, note that 
an axion would feel a CP-violating spin-dependent inter- 
action in presence of matter given by [23| 



V(r) = hg 8 g p - 



1 



1 



8nmc V Xr 



-r/\ 



(36) 



Here, a denotes the neutron spin and ft is a unit vector 
presumably related to the geometry of the macroscopic 
matter configuration. Integrating this potential over the 
geometry of our mirror-absorber-s yste m gives a Yukawa- 
potential contribution of the form [2J| 



A(j)(z) = -a a -^e-^ + a a 



ft P2X r h -z)/X 

8m 3 



-27ra e ff. • A 2 • G 4 • (pi e- z / x - p 2 e 



-(h-z)/\ 



(37) 
) , 



where we have used that 



a c s. 



16ttG 4 • 



4 • m u 



9s9 P 

he 



(38) 



Note, how the switch in the sign of the two terms in 
eq. (|3"8")) arises from the fact that a ■ n/\a ■ ri\ = +1 for 
n the unit normal on the mirror but a ■ n/\a ■ n\ = — 1 for 
n the unit normal on the absorber. Plugging now eq. Ij38|) 
into eq. ()32|) we arrive at a limit for the dimensionless 
axion coupling strength a a given by 



< 



4 m g 
h 2 p 



1 



A/ t (">-((z) -{z)i) 
SAh 



1 



(39) 



where we used that pi = P2 = p in our case. For X^h this 
leads at 90% confidence level to a A-independent upper 
limit on the axion interaction strength 



ME < 2- 10" 

he 



(40) 
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VII. BOTTOM MIRROR AND REAL 
ABSORBER - SMALL A 



This situation is similar to that of Sect. VI, except 
for two crucial differences: Firstly, due to A <C h we have 
to use full potential of eq. (j25|) instead of its expansion 
in eq. ([26]) . Secondly, the presence of the absorber at 
h = h n + At, the turning point, yields a relatively large 
contribution to A<f>(h n ) since there the exponential factor 
in the Yukawa potential is e (,l " _ ' l)/A = e~ AT/A . This in- 
serted into eq. (fT7|) and using eq. (|20|) to compute AE n 
yields an equation for 5h given by 

A£« - m-g-Sh n 

- m ■ A(j) Abs .{h n ) =0 



MO) , 5 

with: AE^ = -C^^^-maG iP -^ 

Here it has been made use of the exponentially decaying 
tail of a bound state Airy function for z > h n that sup- 
presses the contribution of the absorber attached Yukawa 



potential to AEn ■ Therefore we neglect the absorber for 
the calculation of the energy shift. Calculating for the first 
two states 8 Ah = Sho — 5h\ and reshuffling of this equa- 
tion to extract a then writes as a limit symmetrical in its 
sign which reads exactly the same as the simple one of 
eq. m (!) 



a < 



1 



9 



R 3 ■ 6 Ah 



R ^1V R 

A~ 5 for small A, A < 5 /im 



7T • A 5 • G 4 



P 



(41) 



The fact that the limit is symmetrical rests on neglecting 
the absorber in the calculation of the energy shift. Were 
this taken into account it would render the limit asym- 
metrical. Thus, this limit is shown for the attractive case 
in Fig. [3 (thick red long-dashed line) for A > 1.2 \im. For 
A < 1.2 /xm this limit violates the general perturbativity 
bound discussed in the next Section. As we expect the 
limit - as before - for the repulsive case to be weaker than 
for the attractive case, it would be completely outside the 
general perturbativity bound for A < ft,, which is why we 
do not display it. 



Ah 



(2) 



E 



,(°) 



m n ■ A(j) 



m "-9^o mn-g-ih^-hW) 
~ E ^O^i 



7T 2 • a 2 ■ A 4 • G\ ■ p 2 



m>0 



2. 3 2. jR 3. v / /l (°) /l (0) 

with: Cm ~ ■ 



dz ■ sm(<p (z)) sm(cf) m (z)) ■ e z/x 



K ° _ 1 f 1 r 

— ■ m~ 1/3 , k ~ 0.71 and <fi m {z) = — ■ / du\ — 
R R Jo v R 



hm - u 



(use eq. pTt around z = ) 



48 l/3 



.g 2 • W 



E 



1 



(1 + 3/4™) 1 / 3 



o [(m + 3/4)V3 . £ + i] 4 (m + 3/4) 2 /3 - (3/4) 2 /3 



(42) 



— : cr(A) c±. C(3) : convergent 

I 



VIII. VALIDITY OF PERTURBATION THEORY 



There remains now the question of the range of va- 
lidity of perturbation theory. The limit of eq. (|2"Tj) for 
small A in principle might lead to values of a and A where 
A(f) < Eo/m does not hold any more at least for z = 0. To 
estimate the validity of perturbation theory we will resort 
to comparing the 2 nd order perturbation theory to the I s * 



order (see e.g. [25|). Then AE { 2) < AE { Q L > will provide us 
with a validity limit of perturbation theory. The 2 nd or- 
der perturbation theory is given as seen in eq. (|42[) . The 
sum ct(A) is convergent and can be evaluated numerically 



via integral approximations. This, in turn, then yields the 
validity range of perturbation theory for small and large 
A, respectively, as: 



a < 



4-^ 2 -48 1 /3. Gi-p ' X s ' ff(A) 



<t7a) - 1/A 4 



A small 



g-R 7 l i 



5 • 10 



13 



A> 10 /xm 

(43) 

which is plotted in Fig. [2] (green dash-dotted line). 
The matrix elements above have been calculated using 
the WKB approximation for the states. It can be shown, 
however, that using the exact Airy functions for the bound 
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states in numerically calculating the matrix elements pro- 
duces results, that agree with the ones derived from the 
above approximate matrix elements to within 10% on the 
whole range of A. One should note that we have been using 
here the Yukawa potential A<j) including its constant pieces 



ble systematical effects present in these micro-mechanical 
experiments compared to the systematics at work in our 



to yield the full contributions to AE^ and AEq 1 ' (which 
for AE^ at large A leads to a behaviour of AE^ ~ A~ 2 
instead of A -1 , see eq. (f22j) and the one before). 



IX. CONCLUSION 



The recent observation of quantum states of ultra- 
cold neutrons in the Earth's gravitational field allows one 
to derive bounds on the strength and range of an addi- 
tional force from the experimental fact, that the crucial 
measured parameters of the ground state and the first 
excited state, their vertical extensions (which in turn are 
related to their energy eigenvalues), have been determined 
to be consistent with Newtonian gravity to within a dif- 
ferential positioning uncertainty of 5 Ah — 2.0 \im at 90% 
confidence level. Such an analysis is interesting from the 
point of view of systematic errors. Since the absence of 
electric charge and the weakness of its magnetic moment 
extremely suppresses the electromagnetic false effects a 
neutron is exposed to, bounds on additional force derived 
from neutron experiments can be given in a systematically 
quite clean way. 

The analytical bounds derived here and shown in 
Fig. [2] form a consistency check for numerical analyses 
like the one of [J]. Following from the turning point shift 
criterion developed here, our limits are valid for A > 1 /im 
because of breakdown of the perturbation expansion for 
Yukawa-like forces of smaller range. Within a range of 
A = 3 ... 10 \im we observe a change in the negative power 
of the A-dependency of the limit from -2 towards -1. This 
invokes the interesting fact that the expectation of a turn- 
ing point shift 5h ~ AE holding for small A fails for longer 
ranges of the interaction. For even larger A the presence of 
the absorber leads to a A-independent upper limit on a of 
about 2 • 10 11 in the attractive case for A > 10 fim, which 
for this case above A ~ 10 mm further decreases like A~ 2 . 

Furthermore, note that the bound at large A > 5 /xm 
translates into a bound on the strength of CP-violating 
pseudoscalar couplings of the axion within the (previously 
experimentally unaccessible) astrophysical axion window 
which is g p g s /fic < 2 ■ 10~ 15 for 5 /urn < A < 500 /xm. 

The analytical limit on a and A for very small A < 
1 itm given by the analysis of @ is beyond the scope of our 
analysis due to its non-perturbative nature. Other experi- 
mental limits on extra forces are derived from mechanical 
experiments and can be found, e.g., in 

[1,0,1©, El- Those 
limits are derived from Casimir-force measurements or 
mechanical pendulum experiments. They are significantly 
better in numbers than the one derived here, however, 
one should stress the completely different nature of possi- 



l.xlO 15 



tt l.xio' 2 




Limits, this paper: Large A & real absorber, attractive 
Limits, this paper: Small A & real absorber, attractive 
Limits, Ibis paper: Large ,1 & real absorber, repulsive 
Limit 'no-new-state', Nesvizh. et al. 
Perturbation theory validity domain, eq. (30) 
PVLAS signal, [26] 




A in um 



FIG. 2: Red thick solid: 

Limit for large A with real absorber, attractive case eq.s J29} 
and (3TJ. 

Red thick long-dash: 

Limit for small A in presence of a real absorber, attractive case 
eq. 14H . 

Red thin wide-dash: 

Limit for large A with real absorber, repulsive case (eq. H29H ). 
Black solid: 

The 'no-new-state' limit of ref. 0]. 
Green dash-dot: 

Validity boundary of perturbation theory, eq. j43} . 
Blue vertical dash: 

PVLAS signal [2fi| . We display this, as an axion interaction would 
generate a Yukawa-like interaction with a ~ g s gp |23ll . and the 
photon-axion coupling (if PVLAS sees an axion and not a scalar) 
is - up to a loop suppression - roughly of the same order as the 
direct axion-nucleon coupling we measure. 



neutron experiment. Casimir-force measurements depend 
strongly on the geometry of the experiment and the theo- 
retical treatment of the Casimir effect, which is a difficult 
task. Therefore, these limits extracted from Casimir-force 
measurements are not as rigorous as those shown in Fig. [2] 
Using an elementary and electrically neutral particle like 
the neutron one can provide exclusion limits on additional 
interactions untouched by the false effects of the mechan- 
ical experiments. 
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